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Abstract 

We present an analytical calculation of the leading three-loop radiative 
correction to the p-parameter in the Standard Model in the large Higgs mass 
limit. This correction, of order g^m^/M^, is opposite in sign to the leading 
two-loop correction of order g^mfj/M^. The two corrections cancel each 
other for a Higgs mass of approximately 480 GeV. The result shows that it is 
extremely unlikely that a strongly interacting Higgs sector could fit the data 
of electroweak precision measurements. 
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1 Introduction 



The Standard Model of electroweak physics is in good agreement with the data. The 
only particle of the Standard Model, that has not been detected so far is the Higgs 
boson. Direct searches give a lower bound of m# = 114.4 GeV fj]. The precision 
of present day experiments even makes it possible to put limits on the Higgs mass 
through its influence in radiative corrections. At the one-loop level, such radiative 
corrections grow logarithmically with the Higgs boson mass [2]. Fits to the data 
imply a light Higgs boson. However the situation is not completely satisfactory, as 
there is difference of about 2.9 a between the two most precise determinations of the 

effective electroweak mixing angle, sin 2 i-e. the one based on the measurement 

of the 6-quark forward-backward asymmetry ^4pg at LEP on one hand, and the 
one based on measurements of the leptonic asymmetry parameter Ai at SLD on 
the other [3J Hj. The value of the Higgs mass preferred by the b-quark data is 
around 0.5 TeV, while the leptonic asymmetry data and the H^-boson mass point 
to a value which is slightly below the lower bound from the direct searches (5J E] • 
With the recent measurements of the W-mass and top mass from the Tevatron this 
is well within statistics. Still there is the logical possibility that the Higgs boson 
is very heavy (« 1 TeV) and strongly interacting. Since the Higgs self-coupling 
grows like m 2 H higher loop effects can play a role and cancel against the one-loop 
leading log(m#) effects. At the two-loop level such radiative corrections have been 
calculated, also allowing for anomalous Higgs-boson self-couplings [HE]. Inclusion 
of these two-loop corrections does not make it possible to fit the data with a heavy 
Higgs boson. However in a recent paper it was shown that the two-loop correction 
is accidentally anomalously small and therefore important effects might first appear 
only at the three-loop level j^j. The situation can only be clarified by an explicit 
three-loop calculation. This three-loop calculation for one of the precision variables, 
the so-called p-parameter, is the subject of this paper. 

The electroweak p-parameter is a measure of the relative strengths of neutral 
and charged-current interactions in four-fermion processes at zero momentum trans- 
fer pUj- After defining the Fermi constant Gp by means of the effective interaction 
that describes muon decay, 

£ CC = ~ Kl" (1 + 75) /i] [e 7/1 (1 + 75) Ve\ (1) 

one can define p and the sine of the weak mixing angle sw = sin 8w as the parameters 
that appear in the effective interaction that describes the scattering of neutrinos and 
anti-neutrinos by electrons, 

C NC = -0 [lyy" (1 + 75) v,\ [S7, (1 - 44 + 75 ) e] . (2) 
In the Standard Model, at tree level, it is related to the W and Z boson masses by 



where cyy = cosOw This relation gets modified by radiative corrections 



which are sensitive to the existence of heavy particles in the Standard Model, in 
particular the top quark and the Higgs boson. The dominant contributions in the 
large top mass limit come from corrections to the W and Z boson propagators 
involving t and 6-quark loops. They are proportional to m 2L , where L is the number 
of loops, and have recently been calculated up to three loops jTTl [TJ] . 

Another limit one can consider is the large Higgs mass limit. Here, one might 
expect to find corrections of order m 2 ^, since the W and Z self-energies contain 
terms of this order, but it turns out that such terms are screened in low energy 
observables such as the p-parameter [2]. The leading one- loop terms in Ap depend 
logarithmically on m# [T3|, and the leading two-loop correction is proportional to 
m 2 H [TIE]. In this paper, we present the leading three-loop correction Ap®, which 
is of order m A H . 

Radiative corrections to four-fermion processes include self-energy corrections to 
the gauge boson propagators and vertex and box corrections. In principle, all of 
these can affect p and should be taken into account. However, if the fermion masses 
are neglected, the Higgs boson does not couple directly to the fermions. In this case, 
the one-loop vertices and boxes do not depend on m H . Two and three-loop vertices 
and boxes can only contribute to the leading order term of Ap® if they contain 
one-loop subgraphs of order m 2 H or two- loop subgraphs of order m 4 H . However, it is 
possible to choose a renormalization scheme in which all terms of order m 2 H (mjj) 
in the relevant one-loop (two-loop) subgraphs are absorbed into the renormalization 
factors [7J [T3]. In such a scheme, vertices and boxes do not contribute to Ap at 
the leading order in m#, and one is left with just the corrections coming from the 
transversal W and Z self-energies, 

SP(0) S^(0) 

In section |21 we describe the calculation of the bare three-loop gauge boson 
self-energies, which we calculate in the full electroweak Standard Model. The cal- 
culation is simplified by the fact that they are only required at external momentum 
p = 0, which reduces the diagrams involved to vacuum diagrams. Diagrams con- 
taining fermion loops are omitted, since after renormalization, they do not give any 
contributions of order mjj. We simplify the calculation further by systematically 
expanding all diagrams in powers of m# for m# 3> My/,Mz, using the method 
of asymptotic expansions |16j . In this way, all terms are factorized into integrals 
depending only on the large scale ttih on the one hand, times integrals that depend 
on the small mass scales My/ and Mz on the other. The factorized expressions are 
then reduced to a set of independent master integrals using integration by parts 
identities. This reduction is performed exactly in d = 4 — e dimensions. Higgs tad- 
pole contributions and would-be Goldstone boson self-energies, used to determine 
the renormalization constants, are calculated by the same method. An attractive 
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feature of this approach is that it enables one to check Ward identities and observe 
cancellations explicitly at the level of master integrals. 

The renormalization is discussed in section 01 and the final result is presented in 
section HJ 



with 



2 Large mass expansion and reduction to master 
integrals 

All our diagrams are generated by the program QGRAF [T7j . The rest of the calcu- 
lation is done mainly in FORM |18j . For the three- loop one particle irreducible W 
and Z self-energies, for example, there are 104340 and 82985 diagrams, respectively. 
They can be divided into 80 different topologies. After setting the external momen- 
tum p — 0, the diagrams can be expressed in terms of scalar vacuum integrals. We 
use the following notations 

h{m\, m 2 2 , m 2 3 , m\, m\, m\ ; n x , n 2 , n 3 , n 4 , n 5 , n 6 ) = 

J d%d d k 2 d% P(h; mi ) ni P(k 2 ; m 2 ) n2 P(k 3 ; m 3 ) n3 

P(h + k 2 ; m 4 ) n4 P(k 2 + k 3 ; m 5 ) n5 P(h + k 2 + k 3 ; m 6 ) ne (6) 
h{m\, ml, ml, ml, m 2 5 ; m, n 2 , n 3 , n 4 , n 5 ) = 

d d k 1 d d k 2 d d k 3 P(h; mi ) ni P(k 2 ; m 2 ) n2 P{k 3 ; m 3 ) n ' A 

P{k x + k 2 ; m,rP(k 2 + k 3 ; m 5 ) n5 (7) 
h{m\, m 2 2 , m 2 3 , m\;n x , n 2 , n 3 , n 4 ) = 

J d%d d k 2 d d k 3 P(h; mi ) ni P(k 2 ; m 2 ) n2 P(k 3 ; m 3 ) n3 

P(h + k 2 + k 3 ; m 4 ) n4 (8) 
h{m\, m 2 2 , m 2 3 ; n 1 , n 2 , n 3 ) = 

J d% d% P(k x ; mi ) ni P(k 2 ; m 2 ) n2 P(k x + k 2 ; m 3 ) ns (9) 
h(m\ ; nx) = 

J d d h P(k t ; mi ) ni (10) 

1 



P{k ; m) 



k 2 + m 2 

3 



These integrals correspond to the diagrams presented in Figure H In addition to 
the integrals (J5j)- (fTUj) . integrals with scalar products ki ■ kj in the numerator appear. 

In general, the propagators in the diagrams can depend on three different non- 
zero masses: m#, Mw and Mz- The Higgs mass also appears in the scalar 3 and 
4-point vertices via A = g 2 m 2 H / \AM W ). We extract the leading m# dependent terms 
by performing an asymptotic large mass expansion, considering tuh to be large and 
Mw and Mz to be small. Here, we describe the procedure using the language of 
expansion by regions [TEj . 

The expansion is constructed by considering different regions in loop momentum 
space, distinguished by the set of propagator momenta which are large or small in 
those regions. In each region, a Taylor expansion of all propagators in the small 
masses and in the small momenta of that region is performed. Typically, the expan- 
sion generates extra scalar products of loop momenta in the numerator and higher 
powers of denominators, as compared to the original diagrams. However, in each 
term, the dependence of the integrand on the large and small parameters is factor- 
ized. The resulting expression is then integrated over the whole loop momentum 
space. 

For the three-loop vacuum topology Iq, there are 15 regions in loop momentum 
space to consider. 

1. The region where the momenta in all propagators are large. In this case, the 
Taylor expansion yields three-loop vacuum integrals in which all masses are either 
zero, or equal to m H . 

2. Six regions where one momentum, e.g. k±, is small, while the others are large. 
Here, the Taylor expansion leads to products of two-loop vacuum integrals depend- 
ing on run, times one-loop vacuum integrals depending on Mw or Mz- 

3. Three regions where the momenta in two non-adjacent propagators, e.g. k\ and 
ks, are small, and the momenta in the other propagators are large. 

4. Four regions where the momenta in three propagators that are connected to a 
common vertex, e.g. k\, k% and k\ + k 2 , are small, and the momenta in the other 
propagators are large. In these regions and in the regions of type 3, the Taylor 
expansion leads to products of one- loop vacuum integrals depending on m H , times 
two-loop vacuum integrals depending on Mw and Mz- 

5. The region where all momenta are small. In this region, the Taylor expansion 
gives three-loop vacuum integrals depending only on Mw and Mz- These two-scale 
integrals are the most complicated ones that appear in the expansion. So far, only a 
few of the corresponding master integrals have been calculated in the literature • 

Not all regions give non-vanishing contributions for all diagrams. In many cases, 
after the Taylor expansion we are left with scale-less integrals, which are zero in 
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Figure 1: The three-loop topologies J 6 , J 5 and J 4 , and the two-loop topology I 2 

dimensional regularization. Which regions actually do contribute depends on the 
distribution of large masses in the diagrams concerned. For example, in case 2 above, 
the region where k% is small only gives a non-vanishing contribution when mi is a 
small mass. 

Obviously, the factorizable integrals coming from regions 2, 3 and 4 are easy to 
deal with, since they are simply products of one and two-loop vacuum integrals. The 
Iq integrals coming from region 1 can be classified into ten different kinds, according 
to the distribution of massless and massive denominators they contain. For two of 
these categories, we follow the integration-by-parts j2Hl reduction algorithm obtained 
by Broadhurst [21], to reduce them to master integrals. Reductions for the eight 
other categories have been indicated by Avdeev [221, but here, we prefer to use our 
own reduction routines. As a cross-check, we have also performed the reduction of 
these single-scale integrals using the automatic integral reduction package AIR j2S] ■ 
The master integrals themselves are known, some in terms of V functions, the others 
as expansions in e |2H I2S] ■ 

The three-loop integrals from region 5 appear to present a more difficult problem, 
since some of them depend on two mass scales: M w and M z . It turns out that, 
in the three- loop gauge boson self-energies up to the order m A Hl they are all of the 
I4 topology (sometimes with scalar products in the numerator). However, by using 
symmetry relations between such integrals, we find that they cancel out of the gauge 
boson self-energies at this order in m#, provided we sum over all diagrams. 

The gauge boson self-energies we need are the "full" self-energies, consisting of a 
one-particle irreducible part plus terms containing Higgs tadpole insertions, which, 
in turn, include tadpole insertions themselves, as shown in Figures 121 and El 

We have checked several Ward identities that are satisfied by the unrenormalized 
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Figure 2: The full WW three-loop self-energy, the light-grey blobs represent the 
sum of one particle irreducible one, two and three-loop diagrams, the dark-grey 
blobs represent the full one, two and three loop tadpoles. 

self-energies. In particular, we have checked that the photon self-energy at zero 
momentum Y J ^ A '' i ~ loop ^ uU : (0) vanishes at order rrijj. For the would-be Goldstone 
bosons we have checked that £^.3-i«>p,/««(o) and £^V ,3-j oop ,/ u M(q) both vanish at 
order m 6 H . These identities are only valid for the self-energies including tadpole 
contributions. 

3 Renormalization 

We write the Standard Model Lagrangian without fermions as C = Ci nv +Cf ix +Cpp. 
The invariant part is 



r. = --W a W 



1 



B^B^ - (D^ (£>"$) - -A($ t $) 2 - fi $t$ (11) 
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Figure 3: The full two loop Higgs tadpole. 

where W£ v and B^ v denote the curvatures of the SU{2) L and U(l)y gauge fields 
W° and B», 

is the Higgs doublet, and 

= U - Y w ^ a - 4^) $ (i3) 

its covariant derivative. The fields ± , W^, and are defined by 

± = -^{^Titf) , (14) 

W? = ±=(WlTiWl) , (15) 

= cwWl-swB^ (16) 
= swWl + cwB^, (17) 
and the coupling constants g' and g are related to each other by 

9 = ~9 — • (18) 

Cw 



The gauge fixing term is given by 

^ = -c + c--Uc z y-Uc A y no) 



with 



C± = -dpW** + ZMwjt (20) 

C z = -d^ + £M z 4>° (21) 

C A = -d^A" (22) 

In equations (J2U|) - (J2T|) . we have introduced a gauge parameter £. The standard 
't Hooft Feynman gauge fixing term corresponds to £ = 1. Finally, the Faddeev 
Popov Lagrangian C FP is derived from the variation of Cf ix under gauge transfor- 
mations. 

Because the tree level p-parameter is independent of the parameters in the La- 
grangian, Ap® is not affected by any three-loop counterterms. Therefore, we only 
need to renormalize the model up to two loops. Our renormalization scheme is 
similar to the one used in ref. [7]. However, in line with our calculation of the un- 
renormalized diagrams, the renormalization factors are systematically expanded in 
powers of m#, and only terms that can give a contribution of order m A H to A// 3 ) are 
retained, i.e. terms of order m 2 H in one-loop counterterms, and terms of order mjj 
in two-loop counterterms. This implies that we do not renormalize g or cw at all. 

In the Lagrangian, we express A, fi, v and Mz in terms of the four independent 
parameters g, cw, niii and M\y by ji = — \rn 2 H) X = g 2 m 2 H / '(4M^), v = a/2^ 41 , 
and M z = M w /cw Then, the masses m H and M w , and the Higgs and would-be 
Goldstone fields are replaced with bare masses and fields 

(23) 
(24) 
(25) 
(26) 
(27) 

The gauge parameter £ is renormalized in such a way as to compensate the effect of 
the above renormalizations in the gauge fixing term: 

i^ = Z H x ZlH (28) 

We choose the renormalized £ to be equal to one. Each Z-factor is written as 

Z = \- - <5 (2) . (29) 

The renormalization constants Zm and Zh are fixed by imposing conditions on 
the renormalized W and (f) self-energies 

h T \ p 2 =0 ~ (J [SO) 



(0) 

m H 


7 


m H 


M (0) 


= Zm 


M w 


H (0) 


= z H 


H 


£±,(0) 


= Z H < 




00,(0) 


= Zh< 
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and 



_d_ 

dp 2 



,ren I 



lp 2 =0 







(31) 



where the notation X ~ means that X does not contain any terms of order m 2 H at 
one loop, or of order mjj at two loops. These renormalizations remove all the m 2 H and 
mjj terms from the one and two-loop gauge boson self-energies, the self-energies, 
and the mixings between 0's and gauge bosons. The Higgs mass renormalization 
constant Z m is fixed by demanding that 

ReZ HH > ren \ 

. (32) 



m 



H 



Solving equations (j3lH) - (j3~Tj) . we find the following expressions in terms of vacuum 
integrals. 

^2 



r(l) 



r(l) 
'M 



r(2) 



g 



h(m 2 H ;l) 



h(m 2 H ;l) 



(e-4) 

6-e 
4(e-4) 



(33) 
(34) 



-l 2 



z(27r) d M^ 



13 3e 
+ 



+ 



+ m H I 2 (m H , m H , m H ; 1,1,1)— + — + 



128 64 16 (e- 4) 8 (e - 4) 



3s 9 



64 64 8 (e - 4) 



+ m^J 2 (m^, 0,0; 1,1,1) (— + — + 



11 



+ 



1 9 r ( 2 

I it* 

NC? A#2 



;(27r) d M, 



J^m^; !)(- + § + 



IV 



64 64 4 (e - 4) 
1 



4 8 (e - 4) 



(35) 



r(2) 
; A1 



h(m 2 H ;l) 



13 3e 1 
32 ~ 32 + 8 (e - 4) 2 



+ I 2 (m#, m#, m#; 1, 1,1) - 



3e 3 
32 + 8 



+ m 2 H I 2 {m 2 H , 0,0; 1,1,1) -- + - + 



+ 



(e-4) 



+ 



1 9 T ( 2 

NC? A#2 



z(2 7 r) d M 1 



I^m^; 1) l- 7 + 



iv 



32 16 4(e-4) 
1 



4 



(e-4) 



(36) 



(2) 

In order to determine the renormalization constant 5 m , we use analytical results 
for the two-loop on-shell Higgs self-energy from ref . [2H| . Some care is needed here, 
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because the bare parameters M wo and m m used in ref. j2H] do not correspond 
exactly to our bare parameters M$ . This is due to the fact that we do not 
introduce a counterterm 8v 2 for the Higgs tadpole, but instead, explicitly include 
tadpole contributions. Their rriHo/ m H corresponds to our Z m Z H /Z M . Taking this 
difference into account, we find 

^(-1^81088+ IWS+l^ + f^C-g} , (37) 

5 m = + 



16tt 2 / V 4vr J ' \ ' 2J M$y \ 32 e 2 ' e V 64 " ' ~ 64 
SOT 27 ^, 57 n 39 . 261 , 9 ^ . 63 

-7TA 



256 128 



rV31og3 + —-7i\/~3 — — 7r C — ^vr 2 - 1^3 67+^(3)1 
6 32 32 512 32 32 sv ' J 



(3f 



The correction to the p-parameter can now be written as 

v-iZZ,3— loop,ren (r\\ yiWW,3— loop,ren /<->\ 

ApOO = Lr . W - ^ -5 ^ . (39) 



l z lv± w 



The counterterms 8^\ 8^ 2 ' are substituted into equation (|3*9*|) in two steps. In the 
first step, 8jj, 8g , 8j$ and 8$ are replaced with the expressions lJ33|l — (pffijl . (keeping 
8m and 8m as symbols). This step is done before performing any expansions in e. 
At this point, all master integrals depending on My/ or Mz, which originate from 
small momenta regions, cancel out exactly, so that only master integrals depending 
on m H are left. In the second step, the counterterms 8m and 8m are replaced with 
(|37J) -()38J ) and the master integrals are expanded in e, using results from refs. [2"%| l2"5]. 
In the expansions, singular terms of order j = 1, 2, 3, 4, appear and cancel each 
other. The order e term of 8m is not needed, because the coefficient of 8m in 
is finite in the limit e — > 0. Similarly, 8m is not needed beyond the term of order e 
given in eq. (pT7j) . 

As a check on the renormalization procedure, we have verified that the renor- 
malized three-loop photon self-energy, the photon- Z mixing self-energy, and the 
self-energies vanish at zero external momentum. 

4 Results and conclusion 

Finally, combined with the one-loop ^3] and two-loop [7] terms, the three-loop 
correction to the p-parameter reads 

p = 1 + Ap« + Ap^ + Ap® , (40) 
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with 



3 g 2 sly , / in 



4 16tt 2 cl 



log 



(41) 



9 Z 



16vr 2 

2 



Ap 



(3) 



167T 2 
2 



16vr 2 



2 2 2 
g VK m H 

r 2 M 2 
C W m W 

2 2 2 
r 2 /If 2 

3 2 4 
r 2 /If 4 



21 9 
64 ' 32 



+ T^^VS + -r-TV 



3 2 ?CV3 



32 



0.1499) , 
21 729 



(42) 



R 3391 2 9 

7TV3 7T 7TO 

512 512 4608 16 



9109 4 99 R 

7r H v 3 log 3 C 

69120 16 S 



1577 
2304 

— ^Ls 3 (2 7 r/3)-— + — C(3) 



167T 2 



3 2 



Mi 



32 

1.7282) . 



(43) 



The constants appearing in Ap^ are defined by 



U 3 ,i 



v 3ll 



^C(4) + ic(2)log 2 2-llog 4 2-Li 4 (\ 



0.11787599965 

-l) m cos(27m/3) 



(44) 



E 



m 3 n 



m>n>0 

-0.03901272636 



C 



Ch (V 3 ) 



The log-sine integral is defined by 

Ls 3 (9) = 





log 2 


2 sin — 


/o 




2 



(45) 
(46) 

(47) 



Some numerical values are shown in Table 1 and in Figure BJ where we have used 

e 2 47ra 

72~ = T2 - 
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(4£ 



mil /Mw 


Ap« 


r 


r 


2 


- 0.00078 


1.14 10" 6 


-1.33 10" 7 


5 


-0.0018 


7.14 10" 6 


-5.20 10" 6 


6 


-0.0020 


0.000010 


-0.000011 


7 


-0.0022 


0.000014 


-0.000020 


8 


-0.0024 


0.000018 


-0.000034 


9 


-0.0025 


0.000023 


-0.000055 


10 


-0.0026 


0.000029 


-0.000083 


15 


-0.0031 


0.000064 


-0.00042 


20 


-0.0034 


0.00011 


-0.0013 


25 


-0.0036 


0.00018 


-0.0032 


26 


-0.0037 


0.00019 


-0.0038 


27 


-0.0037 


0.00021 


-0.0044 


28 


-0.0038 


0.00022 


-0.0051 


29 


-0.0038 


0.00024 


-0.0059 


30 


-0.0038 


0.00026 


-0.0067 



Table 1: Corrections to p as a function of mn/Mw 

for the weak coupling constant, with a = 1/137 and = 0.23. While Ap( 3 ) is very 
small for low values of m#, it soon becomes more important than Ap( 2 \ The two 
contributions cancel each other at m# « 6Mw « 480 GeV. Ap( 3 ) becomes equal to 
ApW for m H ^2 TeV. 

Comparing our result for the large Higgs mass limit with the ones of the large 
top quark mass limit JT] , we can say that for m# « 500 GeV, Ap® is already larger 
than the three- loop pure electroweak correction of order M 4 6 , which only grows as 
m 2 H log(mj I /Mt ) in the large m# limit. However, it is still very small compared 
with the mixed electroweak/QCD correction term of order a s M^. 

The original question that motivated this calculation was, whether inclusion of 
the three-loop corrections with strong interactions could lead to an effect mimicking 
the one-loop effects of a light Higgs boson. The result of this calculation shows 
that this is highly unlikely. The sign of the three-loop correction is the same as the 
one-loop correction. Therefore with increasing Higgs-mass the three-loop term only 
makes the effect grow faster, instead of the three-loop term partially cancelling the 
one-loop correction. Therefore the presence of a strongly interacting heavy Higgs- 
sector appears to be extremely unlikely, and the data can indeed be used as a strong 
indication for a light Higgs boson sector. 
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Figure 4: The combined effect of Ap^, Ap( 2 ) and A// 3 ) as a function of ttih/Mw. 
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